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We develop an effective multiorbital mean-field t-J Hamiltonian with realistic tight-binding and exchange pa- 
rameters to describe the electronic and magnetic structures of iron-selenide based superconductors AxFQ2-ySQ2 
for iron vacancy doping in the range < y < 0.4. The Fermi surface topology extracted from the spectral 
function of angle-resolved photoemission spectroscopy (ARPES) experiments is adequately accounted for by 
a tight-binding lattice model with random vacancy disorder. Since introducing iron vacancies breaks the lat- 
tice periodicity of the stochiometric compound, it greatly affects the electronic band structure. With changing 
vacancy concentration the electronic band structure evolves, leading to a reconstruction of the Fermi surface 
topology. For intermediate doping levels, the realized stable electronic structure is a compromise between the 
solutions for the perfect lattice with y = and the vacancy stripe-ordered lattice with y = 0.4, which results 
in a competition between vacancy random disorder and vacancy stripe order. A multiorbital hopping model is 
parameterized by fitting Fermi surface topologies to ARPES experiments, from which we construct a mean-field 
t-J lattice model to study the paramagnetic and antiferromagnetic (AFM) phases of K0.8Fe1.6Se2. In the AFM 
phase the calculated spin magnetization of the t-J model leads to a checker-board block- spin structure in good 
agreement with neutron scattering experiments and ab-initio calculations. 

PACS numbers: 78.70.Dm, 71.10.Fd, 71.10.-w, 71.15.Qe 



I. INTRODUCTION 

The discovery of high-Tc superconductivity in iron- 
selenide superconductors Ax¥Q2-y^^2 (A =K, Rb, Cs, Tl, 
Tl/K, Tl/Rb)^ has generated increased excitement in the 
search for high-temperature superconductivity. It provides a 
new opportunity to understand the underlying physics of iron- 
based superconductors. The new family of compounds has 
several unique features: (i) Superconductivity (for x ^ 1 
and y ~ 0.12 — 0.3) emerges in proximity to an insulating 
phase^^ (for x ^ 0.8 and y > 0.4), instead of a poor metal 
as in other iron-based parent compounds. For the iron defi- 
cient compounds with y > 0.4, there is mounting evidence for 
the existence of iron vacancy ordered superstructures stabi- 
lized with a stripe-like antiferromagetic (AFM) state.'^ This 
raises the interest in the possibility of the insulating phase be- 
ing driven by the Mott localizatiorP^EH due to the reduction 
in kinetic energy .'^ In particular, the compounds with x ^ 0.8 
and ~ 0.4 are of special interest for the formation of a pecu- 
liar vacancy order (so called a/5 x >/5 superstructure) as well 
as a block- spin antiferromagnetic (BAFM) state. ^^ ^^ (ii) The 
end members of the series AFe2Se2 (x = 1 and ^ = 0) are 
heavily electron doped (0.5 electron/Fe) relative to other iron- 
based superconductors (such as LaOFeAs, BaFe2Se2, FeSe 
etc.). Band structure calculationi^^^ for these end com- 
pounds show only electron pockets that are primarily located 
around the M point of the Brillouin zone (BZ) as defined 
for a simple tetragonal structure. Indeed a series of angle- 
resolved photoemission spectroscopy (ARPES) experiments 
has been performed on r ecent ly discovered superconducting 
compounds Ax¥Q2-y^^2^^^ ^ common feature is the pres- 
ence of electron pockets around the M point in the Brillouin 
zone and a marked absence or near absence of a hole pocket 
at the F point. With regard to electron pockets, the FeSe-122 
family is similar to the isostructural FeAs-122 family, while 



they differ with respect to the hole pocket, which is consid- 
ered to be essential for certain interband pairing models of 
superconductivity. 

These unique features raise the hope to gain new insights 
into the mechanism of iron-based superconductivity by study- 
ing the FeSe-122 family. However, special care must be taken 
in the interpretation of these results due to the complicated 
real-space structure of highly iron deficient compounds. The 
real-space structure for different Fe compositions (0 < y < 
0.4) is quite intricate and resembles more that of an alloy 
than of a lightly doped crystal. For example, although both 
74Fe2Se2 (y = 0) and 74o.8Fei.6Se2 (y = 0.4) compounds 
have perfect lattice periodicity (albeit the latter one exhibits 
striped vacancy order) in the iron layer, the lattice structure for 
compounds with < y < 0.4 can be thought of as a superpo- 
sition of both lattices, which forms either a random- vacancy 
lattice or phase- separated lattice with vacancy stripe order. 
More generally, the serious problem of iron deficiency is that 
the introduced random-disorder scattering centers in the iron 
layer destroy the translational periodicity, thus rendering the 
wave vectors of the Bloch wave functions as 'bad' quantum 
numbers to describe electron motion. Therefore, when inter- 
preting the electronic structure as measured by ARPES, which 
probes the momentum space, a real-space electronic structure 
approach must be developed to account for the strong disor- 
der. 

In this paper, we present a systematic study of the evo- 
lution of the normal-state electronic structure and magnetic 
properties with vacancy doping based on a real- space descrip- 
tion. The Fermi surface topology extracted from the spectral 
function of the ARPES experiments is adequately accounted 
for by a tight-binding lattice model with random vacancy or- 
der at the Fe sites. We find that introducing Fe vacancies 
breaks the lattice periodicity and drastically affects the elec- 
tronic band structure and Fermi surface topology. The evolu- 
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FIG. 1. (Color online) Different cuts (unit cells) of the real-space 
lattice of FeSe-122 (a) and corresponding unit cells (same color map) 
in the Brillouin zone (b). The arrows indicate the block- spin AFM 
structure. For unit cells with cuts "-b-" and "-C-" the Se atoms in 
the layers above and below are shown. In the vacancy-ordered state 
(y = 0.4) vacancies phase- separate to form stripes. 



tion in the band dispersion results in a noticeable reconstruc- 
tion of the Fermi surface (the technical detail is provided in 
the Appendix). As a consequence, for intermediate iron va- 
cancy concentrations, the realized stable electronic structure 
is a compromise between the solutions for = (perfect lat- 
tice) and y = OA (stripe ordered lattice), resulting in a com- 
petition between vacancy random disorder and vacancy stripe 
order. Based on such a parameterized hopping model, the con- 
structed mean-field t-J lattice model gives rise to a checker- 
board block- spin structure for K0.8Fe1.6Se2, which is in good 
agreement with neutron scattering experiments and ab-initio 
calculations. 

The outline of this paper is as follows. In Sec.[ll|we for- 
mulate a tight-binding t-J model Hamiltonian and introduce 
within the mean-field approach the Bogoliubov-de Gennes 
(BdG) equations. In Sec. |lll| we discuss the Bloch wave 
function formulation of multiorbital electron hopping to ob- 
tain a single set of model parameters for the kinetic energy 
part of the Hamiltonian by fitting the electronic structure and 
Fermi surfaces of Ka^Fe2Se2 with a perfect lattice structure. 
In Sec.|lV]we discuss the electronic structure of a random va- 
cancy lattice by introducing an auxiliary impurity scattering 
approach in the unitarity limit. In the case of the supercell 
calculations, the results are in good agreement with the Bloch 
wave function method. In Sec.|V]we present our calculations 
of the magnetic structure, which agree well with the neutron 



scattering measurements. The summary is given in Sec. VI 
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FIG. 2. (Color online) Schematic pictures of the effective t-J Hamil- 
tonian. (a) The effective tight-binding model, where ti is the nearest- 
neighbor inn) hopping, ti (ts) are the next-nearest-neighbor innn) 
intra-orbital hopping terms due to up (down) Se atoms and t^ is the 
nnn interorbital hopping, (b) The effective exchange interactions, 
where J\ and J2 are the nn and nnn interactions inside each block, 
while J'x and J2 are the nn and nnn interactions between blocks. 



mediation of upper and lower As atoms could lead to differ- 
ent hopping terms between iron atoms in the iron layer. It 
has been noted that this Hamiltonian has a lower S4 point 
group symmetry with respect to the As atoms as compared to 
the tetragonal C4 symmetry of the bulk crystal, which should 
give rise to additional exotic magnetic and superconducting 
phases.!^ Since this tight-binding model was introduced, sev- 
eral successful studies have be en pe rformed^^ to describe 
ARPES,^^^^ magnetic structures, I^^^H phase diagrams,^ and 
vortex core and spin susceptibility in RPA calculations. ^^^1^ 

Figure [T] shows the schematics of various lattice config- 
urations with real- space unit cells used in this work for 
Aa;Fe2-2ySe2. Recognizing the crystallographic similarities 
between the isostructural FeSe-122 and FeAs-122 compounds 
and that the electron pocket at the M point is not only a com- 
mon but also main feature in the heavily electron-doped re- 
gion, < y < 0.4, we use the same tight-binding model for 
the electron hopping as in Refs. [25land[33] We account for 
the reported electronic band structure of the perfect lattice of 
Ka;Fe2Se2,^^ shown in Fig.jSJa), by proposing a modified set 
of hopping parameters (^1,^2,^3,^4) = (1,1,-2,0.08). 

In our approach, we use for simplicity the same hop- 
ping parameters for K0.8Fe1.6Se2 as for all other Ka;Fe2Se2 
compounds. The unit cell of K0.8Fe1.6Se2 is modified to a 
>/T0 X >/T0 area, see Fig.[TJa) with cut "-c-" for the unit cell, 
due to the periodic vacancy order along stripes at the doping 
concentration y = 0.4. In this real-space unit cell there are 
a total of 10 sites, namely, 8 iron atoms and 2 vacancies. We 
start with an effective lattice model for K0.8Fe1.6Se2 by in- 
cluding the hopping and exchange interaction H'^ terms: 



II. MODEL AND FORMALISM 

In this work, we adopt the tight-binding model Hamiltonian 
that successfully describes the structurally related FeAs-122 
superconductors. Zhang and coworkers^^ suggested that the 
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and 



The hopping parameters on the lattice are defined as 
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where i (j) are site indices, and fi {u) are orbital indices cor- 
responding to dxy or dyz wave function orbitals and to is the 
chemical potential. The expressions < ij > (<C ij ^) 
and < ij >' ij denote intra- (inter)-block nearest- 
neighbor inn) and next-nearest-neighbor (nnn) hopping pro- 
cesses, whereas t^PP^^ (^^lower^^ indicates the hopping term 
crossing over the upper (lower) Se atom as shown in Fig.|2ja). 
If we further approximate the exchange interaction to be of 
the Ising type then only the component is involved 



(5) 



and the interaction term can be expressed in mean-field ap- 
proximation by 



(6) 



In Fig.[2[b) the nn intra- (inter)-block exchange term Ji (J() 
and nnn intra- (inter)-block exchange term J2 ( J2) are illus- 
trated. We can now construct the corresponding mean-field 
Bogoliubov-de Gennes (BdG) matrix equation on a lattice: 




^^"^ =E„ 



(7) 

Here the single-particle Hamiltonian iJj^^v o- is expressed by 

5p (8) 

where a correspond to ±1 for spin-up(down) index, S is taken 
with ±x (±y) and ±x ± y, which correspond to the nn and 
nnn real space shift; J (J') is for the intra (inter) block in- 
teraction and the mean-field magnetization per site and per 
orbital is (m^^) = /iB((^i/it) " (^i/ii))- The quasiparticle 
energies En are measured with respect to the chemical po- 
tential. We note that the single-particle Hamiltonian depends 
on spin- and orbital-dependent electron density and pairings. 
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FIG. 3. (Color online) The electronic band structure for compounds 
(a) Ko.8Fe2Se2 (fill factor n = 2.4, to = 1.09 for high electron dop- 
ing) and (b) K0.8Fe1.6Se2 (half-filling n = 2.0, to = 0.005) with 
(ti, t2, ts, t4) = (1, 1, —2, 0.08). The dashed blue line is the loca- 
tion of the chemical potential, (c) Fermi surface of electron pockets 
of K0.8Fe1.6Se2 located near the F point in /c-space corresponding to 
the y/lO X >/T0 unit cell (blue square). a (d) Static spin susceptibility 
for the bare band structure of K0.8Fe1.6Se2. The high intensity spots 
at scattering wave vectors Qv = 7r(|, |) and Qm = 7r(|, |) are in 
agreement with neutron scattering experiments. 



which are given by 

n 



(9) 
(10) 



A.^i. = ^ E«MT^l"i + <.T<" ) tanh ^.(11) 

n 

Here we set the Boltzmann constant ks = I- Finally, the BdG 
equation ^ must be solved self-consistently with Eqns. 
( pT) . Since we consider only normal-state properties in this 
work, we ignore the superconducting pairing term, Fi^^v = 0, 
in the BdG matrix equation. 



III. ELECTRONIC STRUCTURE IN THE 
PARAMAGNETIC STATE 

The most common way to construct the electronic band 
structure of ordered systems is to use the Bloch wave function 
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formulation by In case of the enlarged ylO x ylO unit cell 
with 8 Fe atoms and 2 vacancies, the corresponding k-space 
Hamiltonian is a 16 x 16 matrix, i^fgxie' which is straight- 
forward but tedious to derive. Instead, we propose another 
method based on the impurity problem solution that produces 
exactly the same results and, moreover, provides physical in- 
sight into the breaking of the translational symmetry and how 
it gradually affects the electronic structure with increased va- 
cancy doping. In this approach the vacancy is mapped onto an 
impurity with an adjustable onsite scattering potential Vq that 
varies from (no vacancy) to oo (strong scattering center). 
Also the continuous tunability of the impurity potential pro- 
vides an easy handle on the evolution of the electronic struc- 
ture with scattering strength. Please refer to the Appendix for 
more details on the implementation of the impurity problem 
calculation employed in this work. 

In Fig. |3ja) we show the electronic band structure of 
Ko.8Fe2Se2 for our tight-binding model parameterization us- 
ing the V2 X >/2 unit cell (cut "-b-" in Fig. [l|a)). The cor- 
responding band structure of the vacancy ordered compound 
K0.8Fe1.6Se2 is shown for comparison in Fig.jSjb). The evo- 
lution of the dispersion from the a/2 x a/2 unit cell with four 
bands to the enlarged a/TO x a/To unit cell with 16 bands is 
nontrivial and cannot be obtained by a simple rigid-band shift 
of the chemical potential. Indeed, the Fermi surface topol- 
ogy for K0.8Fe1.6Se2 shows very tiny electron pockets located 
near the F point, see Fig.jSjc), in contrast to Ko.8Fe2Se2. 

In addition to the electronic dispersion, we calculate the 
static spin susceptibility x(q) for the bare band structure of 
itinerant electrons. 



x(q) = X^x/j(q,^^m 0) , 



ij 



X/J = ^(T.[5|(q,r)55(-q,0)]), 



(12) 



(13) 



where N is the total number of unit cells and 6'|(q, r) is the 
Fourier component corresponding to the site and orbital index 
/ in the unit cell. In this compact notation, the orbital wave 
functions have the running super-indices / = (i, /i) and J = 
(j, u). The dynamic susceptibility is given by 

Xij{q,i^m) = ^^Gjc,,/a(K,ia;n) 

Kq! n 

X Gia,Ja (K + q, iuJn + i^m) , (14) 

where the wave vector K is defined in the Brillouin zone cor- 
responding to the a/To X a/To unit cell, ujn = (2n + l)7rT 
and rtrn = 2m7rT are the Matsubara frequencies of the 
fermions and bosons. In standard notation the multiorbital 
lattice Green's functions are given by 



G/;,j;(K,ZCJn) 
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FIG. 4. (Color online) Schematic picture of a 20 x 20 lattice with 
random vacancies. For this configuration the minimum size of a su- 
percell is a 10 x 10 square lattice. This supercell contains the trans- 
lational symmetries of both the ^2 x a/2 and a/5 x a/5 unit cells. 
The random lattice is constructed by randomly removing vacancies 
from the ordered diagonal stripes of phase- separated vacancies. The 
phase separation pattern along the direction (2,1,0) becomes clearer 
for larger superlattices. 



In the calculation of the itinerant spin susceptibility the red 
spots in Fig. [3jd) show high intensity around q = and 
Qa in agreement with neutron scattering experiments.^^ It 
follows from the Stoner criterion that the observed AFM state 
is possibly formed from itinerant electrons of the paramag- 
netic state due to the bare band structure, rather than the ex- 
change interaction of localized spins. Hence it is natural to 
expect that close to the magnetic instability a small driving 
force can break the symmetry of the paramagnetic state and 
induce the long-range AFM state. 



IV. ELECTRONIC STRUCTURE OF RANDOM VACANCY 
LATTICE 

The periodic Bloch wave function formulation becomes 
questionable for lattices with random vacancy disorder in the 
intermediate doping region < y < OA, where a perturbative 
approach to disorder in supercells fails. Here, we describe the 
construction of a random vacancy lattice calculation (RLC). 
First, we construct a 40 x 40 or larger real-space lattice for 
AajFe2Se2 with vacancies phase- separated in stripes along the 
(2,1,0) direction as shown in Fig. |4] Second, we randomly 
add iron atoms, i.e., remove vacancies, on the vacancy ordered 
stripes up to the point that we match the vacancy occupation 
for compounds measured by ARPES. Third, we construct the 
corresponding real- space Hamiltonian with hopping terms 
(^1,^2,^,^4) = (1,1,-2,0.08). Finally, we exactly diago- 
nalize the Hamiltonian and compute its pairs of eigenvalues 
En and eigenvectors (i^^^^, ^i^^) from which we calculate the 
lattice Green's functions. 

When going beyond the dilute limit of random disorder, 
as for example in a heavily doped alloy, the randomness of 
vacancies plays a crucial role for the electronic dispersion. 
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because the rigid-band shift appropriate in the dilute Hmit 
breaks down. Hence the conventional approach of periodic- 
ity and the calculation of the electronic band structure and 
Fermi surface fail. In order to solve this complicated prob- 
lem, we calculate the multiorbital real- space spectral function 
^/j(^) = —^^^Gij{iuJn ^ uj -\- iV), which depends on 
the multiorbital real- space Green's function summed over all 
bands n, 

Gij{iuJn) = G/^, Jt(K = 0, iUn) + G/;,j;(K = 0, iUn) , 

(16) 

where the general form of the spin-dependent Green's func- 
tions on the right-hand side of the above equation has been 
given in Eq. (15). Finally, the Fourier transform of Aij{uj) 
to k- space for the one-iron per unit cell (1-Fe) gives the de- 
sired spectral function, which is the one measured in ARPES 
experiments, 

A{k,uj) = ^A/j(cj)exp{-ik- (Ri -R^)}. (17) 
ij 

Here, i (j) are lattice indices for the location of each Fe atom 
in the 40 x 40 supercell with orbital indices /i (u) and super- 
indices / = (z, /i) and J = (j, u). The double-sum is a short- 
hand notation for = ^ S^jy. Since we considered 
only the bare band structure, we dropped for convenience the 
spin indices in this calculation. 

Next we apply the RLC method to calculate for both dop- 
ing variables x and y the evolution of the spectral function and 
Fermi surface topology for compounds AxFe2-y^^2 • For the 
hopping term, H\ the variables x and y can be mapped onto 
the value of the chemical potential to (x) and vacancy concen- 
tration y. Figure [5] shows the calculated Fermi surface evo- 
lution for changing to{x) and y values on random vacancy 
disorder lattices. These lattices were already large enough 
to self-average the vacancy disorder, thus no further ensem- 
ble average of different random vacancy configurations was 
required. The visualization of the evolution of the spectral 
function helps us to understand when an electron pocket "ap- 
pears" at the F point, namely for y > 0.22, see panels (b)-(d), 
as well as when it disappears for shifted chemical potentials 
t^Q^ = -0.7 and 4"^^ = -1-5 (rows 3 and 4). Our finding dif- 
fers from that discussed in Ref . 34 , where it was claimed that 
the electron pocket at the F point of the FeSe-122 compound 
originates from the BAFM structure. However, our results 
clearly demonstrate the presence of an electron pocket at F in 
the absence of BAFM order and over a large range of vacancy 
concentrations 0.22 < y < 0.4. Indeed Fig. [5] shows that the 
small electron pocket at the F point originates from the same 
regions in the Brillouin zone as for the vacancy stripe-ordered 
compound with Fei.e shown in Fig. |3] These calculations 
demonstrate the power of the RLC method, when randomness 
(doping) strongly affects the Fermi surface. This computa- 
tional method is adequate to understand the evolution of the 
electronic structure for intermediate doping levels in the range 
< ^ < 0.4, where the stable electronic structure is a com- 
promise between the solution for the perfect lattice with ^ = 
and the stripe-ordered lattice with y = 0.4. Alternatively one 
can obtain these results by diagonalizing large supercells with 
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FIG. 5. (Color online) The calculated Fermi surface evolution from 
the Fourier transformed real-space spectral function in the 1-Fe BZ. 
The alphabetical order in each picture represents different vacancy 
concentrations: (a.n) Fei.g, (b.n) Fei.rs, (c.n) Fei.72, (d.n) Fei.7, 
while the numerical index n — 1—4 represents different values of the 
chemical potential. The chemical potential varies from high to low 
as n increases: t[,^\ 4^^) = (0.5,0.33,-0.7,-1.5). The 
calculations in panels (b.l), (c.2) and (d.3) correspond to the com- 
pounds (Tl,K)Fei.78Se2, Tlo.58Rbo.42Fei.72Se2 and K0.8Fe1.7Se2, 
respectively. In panel (a.4) we redraw the different unit cells of 
Fig.[l];b) with 1 X 1 (green), ^2x^2 (red), and ^10 x ^10 (blue). 



TABLE I. The FeSe-122 compounds and their respective vacancy 
doping y, whose Fermi surface topology was fitted to local density 
approximation (LDA) calculations or ARPES measurements. 
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Ko.8Fe2Se2 
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(Tl,K)Fei.78Se2 


ARPES 
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Tlo.58Rbo.42Fei.72Se2 


ARPES 


RLC 
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0.30 


K0.8Fe1.7Se2 


ARPES 


RLC 
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0.40 


K0.8Fe1.6Se2 


LDA 


PLC 
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random lattice Hamiltonians. Of course such a brute force ap- 
proach is time consuming, especially when many calculations 
for different random configurations are required as illustrated 
in the cartoon of Fig.|4] We confirmed numerically our results 
by sampling many random vacancy configurations of large su- 
percells and found that the spectral functions are no different 
from the RLC method presented here for fixed iron concentra- 
tions. In Table |l] we list results for both periodic lattice calcu- 
lations (PLC) and random vacancy lattice calculations (RLC) 
performed for various vacancy concentrations y. 
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FIG. 6. (Color online) (a) Temperature dependence of the magnetiza- 
tion (as defined in the main text) in the BAFM state. The temperature 
T is in units of ti . The inset shows a small part of the lattice used in 
our mean-field calculation (20 x 20 sites with 20% vacancies), (b) 
The zero-temperature LDOS as a function of energy uj in units of ti. 
Shown are the two block sublattices with net spin up (dashed red) 
and down (dashed blue) lines as well as their sum (solid green). A 
gap of 0.4 ti opens at the Fermi level a; = 0. 



V. MAGNETIC STRUCTURE IN THE AFM STATE 

Now that we have fully developed the parameterization 
for the kinetic part of the t-J lattice model for random 
vacancy order with a single set of hopping parameters, we 
can focus on the exchange interaction term H'^ . We start 
with a mean-field calculation for the magnetic state, using a 
set of rescaled exchange parameters extracted from the band 
structure calculations.^^ We rescale the ab-initio exchange 
parameters by a global factor to obtain ( Ji, J2, J{ , J2) = 
(—3.44, —0.36, —1.16, 1.52)^1, which results in better agree- 
ment of the calculated Neel temperature with experiment. 
Within our low-energy effective t-J model, this set of ex- 
change parameters leads to a Neel temperature for the block- 
spin AFM state, Tn ^ 1.4 ti ^ -0.407 Ji, see Fig. |6];a). 
If we assume a hopping integral ti = 35 meV, we obtain 
Tn = 569 K. Indeed this value is very close to the one ob- 
served experimentally (Tn = 559 K).^^ We argue that this 
quantitative result provides further evidence for the domi- 
nance of the kinetic term over the exchange term in the t-J 
Hamiltonian with incipient AFM order due to Fermi surface 
nesting. Furthermore, we calculate the staggered magnetiza- 
tion as defined by 



-y\ 



and its deviation from the standard AFM state may be defined 
as 



(19) 



(18) 



Here Ni^ is the total number of Fe sites and = 
is the spin density on each site where the summation index 
i runs over all lattice sites. In the BAFM state, the system 
breaks the translational symmetry in the v^S x unit cell 
and gives rise to a larger periodicity with a >/T0 x >/T0 unit 
cell. The result for the magnetization as a function of tem- 
perature is shown in Fig. [6ja), which is in qualitative agree- 
ment with the neutron scattering experiments, except for a 
smaller magnetic moment as observed experimentally per iron 
site (3.3/iB)PThe reason for this discrepancy follows from 
our phenomenological two-orbital model, which has only 2 
electrons per Fe site. Therefore, the maximum total moment 
on each site is I^ib- Whereas all 5 electrons of the Fe atom 
seem to participate in the moment formation. Notably our 
result shows an antiferrimagnetic state at finite temperatures 
for T > 0.3ti, where ^ 0, that is, the magnetization 
for the up and down spin blocks are slightly different lead- 
ing to a net ferri magnetization. This difference is related to 
the fact that in the present model an asymmetry for the intra- 
orbital next-nearest-neighbor hopping integrals ^2 vs. t^, has 
been introduced due to the selenium atoms below and above 
the iron layer. We confirmed numerically that for 1^2 1 = 1^3 1 a 
purely antiferromagnetic state exists. It is interesting to note 
that this type of difference in kinetic hopping terms also mani- 
fests itself in the magnetization between the two spin blocks of 
the bipartite sublattices. The existence of antiferrimagnetism 
could be identified by weak satellite peaks at the ferromag- 
netic propagating wave vector Q = in elastic neutron scat- 
tering of K0.8Fe1.6Se2 for temperatures 0.3ti < T < Tn- 
Therefore using the same arguments one should also expect 
an asymmetry in the AFM split satellite peaks of the nuclear 
magnetic resonance (NMR) spectrum. The observation of an- 
tiferrimagnetism could serve as an important evidence for dif- 
ferent Fe 3(i-electron hopping integrals mediated by Se atoms 
buckling on both sides of the Fe layer in the bulk system. 

Finally, we calculate the local density of states (LDOS) 
given by 

Mco) = J2[K^^\H{E„ -co) + K^fd{Er,+co)], (20) 

where the delta- function S{x) is approximated as S{x) ^ 
r/7r(x^ + r^) with the intrinsic quasiparticle broadening pa- 
rameter r. In our calculation, we choose F = 0.002 but note 
that the result does not change qualitatively for small changes 
of F. For the LDOS calculations we use an 80 x 80 sites 
supercell. Figure |6jb) shows the zero-temperature LDOS in 
the BAFM phase on two block sublattices. We note that the 
LDOS on each site of a chosen block is the same. As a conse- 
quence of the exchange interaction a clear gap of size ~ 0.4ti 
opens at the Fermi energy cj = 0, which is in qualitative agree- 
ment with ab-initio band structure calculations for the BAFM 
state. 
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VI. CONCLUSION 

In summary, we systematically addressed within an ef- 
fective two-orbital model the effects of the iron vacancy on 
the normal- state electronic structure and magnetic properties 
in the AxFe2-y^^2 compounds. We determined a suitable 
choice of hopping parameters for Fe Sd-electrons by fitting 
the Fermi surface topology of ARPES data. From the deter- 
mined hopping parameters we calculated the evolution of the 
electronic structure and mapped out the Fermi surface topol- 
ogy for arbitrary Fe vacancy concentration based on a ran- 
dom vacancy disorder model. After the kinetic part of the 
t-J Hamiltonian was parameterized we focused on the special 
case of A0.8Fe1.6Se2, where the Fe atoms form the x ^/5 
vacancy order. Here we studied the magnetic properties by 
including the exchange interactions with relative strengths ex- 
tracted from ab-initio calculations. Our mean-field solution of 
the t-J Hamiltonian reproduced the block- spin antiferromag- 
netic (BAFM) state, in good agreement with neutron scatter- 
ing experiments and ab-initio calculations. Finally, we found 
that the magnitude of the magnetization on the two-block 
bipartite sublattices is at a small variance for temperatures 
O.Sti < T < Tat. This feature is unique to our model with 
different Se-mediated hopping strengths (below and above the 
Fe layer) on the two block spin sublattices, which can be 
tested by future refined neutron scattering and nuclear mag- 
netic resonance experiments in the corresponding temperature 
range. 
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Appendix A: Treatment of the vacancy impurity problem 

We treat the vacancy disorder as an impurity scattering 
problem, since one might anticipate that a vacancy behaves 
like a strong impurity scatterer for a Bloch wave function. 
Therefore, we put an on-site impurity energy on each va- 
cancy site i = iy to mimic vacancy disorder: 



H — Hjiop -\- Himp 



(Al) 



We introduce a new index to describe the Ai(2);^i(2) 
sub-atoms for unit cell cut "-b-" in Fig. [l] tpj = 




FIG. 7. (Color online) The band structure plot in high symmetry axis 
(r ^ X ^ M ^ r) for H = (a), Vo = 2.5 (b), Vo = 100 (c). 



(^Ai' ^A2' ^Li' ^52)^' linearize the fermionic operator, 

CAaiBa),i = (1/ V^) Ek CAa(Ba) ,k X exp(zk • Ri), tO get 

the 4 X 4 hopping matrix in k-space, 
( 



£T,k ^xyM 



\ 





^T,k 
^xy^\<i 
^B,k - M / 



(A2) 



We perform the same linear transformation for the impurity 
term Hi^np in the Hamiltonian. There are four independent 
vacancy ordering vectors as shown in the unit cell cut "-b-" in 
Fig.gfor Qi = 27r(4/5, 2/5), Qn = nQi.n = I • • • 

By construction Hirnp enlarges the basis in the k-space, 

V^k ^ V^L = (V^k,V^k+Qi,V^k+Q2,V^k+Q3,V^k+Q4). In this 
formulation the diagonal scattering matrix V = Iaxa x Vq 
sits in the off-diagonal positions to describe multiple scatter- 
ing between all vacancy ordered states with different Qn, Fi- 
nally, the enlarged Hamiltonian becomes H = Ek V^k ^kV^k 
with 



(A3) 

Note that the impurity potential Vq not only leads to a recon- 
struction of the shape of the band structure, it also creates a 
gap of 2Vb between high-energy bands and the low-energy 
bands. Therefore, we need to shift it back to the origional 
place, Wk^Wl. = Wk^ hox20 x Vq. 

Figure |7] shows the calculated band structure for various 
values of the impurity scattering strength Vq. When Vq = 
0, the band structure exhibits an entanglement of 20 bands, 
which down-folded simplifies to the band structure shown for 
the four bands in the y/2 x \/2 unit cell of Fig. [TJ a). The more 
complicated plot in Fig.jTJa) it is due to the repeated plotting 
of different Qn's in k-space of Mk. However, for finite val- 
ues of Vo a gap separates the upper four upper bands and 16 
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lower bands, shown in Fig.|7jb). Finally, when Vq approaches 
the unitarity limit. Vq oo, as shown in Fig.jTJc), the upper 



four bands are being pushed far above as four independent flat 
bands, while the lower 16 bands form a simpler shape as the 
new periodicity with vacancy stripe order sets in. 
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